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Abstract 

A black brane solution to a Gauss-Bonnet massive gravity is introduced. In the 
context of AdS/CFT correspondence, the viscosity to entropy ratio is found by the 
Green-Kubo formula. The result indicates violation of the well-known KSS bound as 
expected in a higher derivative theory. Setting mass zero gives back the known viscosity 
to entropy ratio dependent on the Gauss-Bonnet coupling, while without Gauss-Bonnet 
term, a nonzero mass parameter doesn’t contribute to the ratio which saturates the 
bound of A. 
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1 Introduction 

For decades, various modifications of the Einstein gravity like Lovelock gravity |1], 
brane world cosmology [2], scalar-tensor theories[3], and the so-called F(R) gravity have 
been proposed to address important problems in cosmology such as the cosmological 
constant, dark energy and dark matter. A common ingredient of these modified theories 
is the massless graviton. So it would be interesting to consider addition of a massive 
graviton to these theories 0IU- To study this new class of gravities it is important to 
find various spacetimes in different setups. One of the interesting setups is inclusion 
of the massive gravity in a higher derivative theory such as the Gauss-Bonnet gravity 
and looking for a black hole solution. This is what we do in section 2. In the context of 
massive Gauss-Bonnet gravity with a cosmological constant, we find a class of solution 
which includes black holes with either spherical or hyperbolic horizons and particularly 
a black brane solution with a flat horizon. 
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On the other hand, it is well known from the AdS/CFT duality [9, TU[ ]] £] that 
quantum gravity, string theory, on an AdSd+ 1 background is dual to a d-dimensional 
CFT, which is a non-gravitating theory. The latter can be described by hydrodynamics 
as an effective theory of QFT at large distances and time-scales [12]. The AdS/CFT 
duality leads to fluid/gravity correspondence in this limit [13], E] 15]. So it is worth 
to study the hydrodynamics of the QFT by the corresponding gravity theory via the 
AdS/CFT duality. The suitable solution to investigate the hydrodynamics is a black 
brane. 

The hydrodynamics equations are laws of conservation of energy and momentum 

mm, 


= 0 , ( 1 ) 

r = (/)+ P )«v+pf. ( 2 ) 

According to dictionary of AdS/CFT, a black brane within the bulk is dual to a fluid 
on the boundary. It also implies that the Einstein equation in the bulk corresponds 
to the conservation of the energy-momentum equation on the boundary and the bulk 
metric is the dual field of the energy-momentum tensor of the boundary theory. 

In the large wavelength limit where the hydrodynamics regime is valid, we expand 
the energy-momentum tensor as follows, 

= {p + p)u^u u + pg^ - a^ (3) 

a ^u = pfiapuf) [rii y aUp + V p Ua ) + ( C _ lr,)g a ^.u] (4) 

P> w = g lw + vPu\ 

where g, a^ and P are shear viscosity, bulk viscosity, shear tensor and pro¬ 
jection operator, respectively [El E] US] ESI [2D]. Here we are interested in the shear 
viscosity which can be derived by Green-Kubo formula [201121]. 

V = Um ^ f dt dx e lujt ([T* ( x ), T* (0)]) = - Jim ^f$G x y * (u, 0). (5) 

Notice that the Green-Kubo prescription is independent of details of the bulk theory. 
Furthermore, the Gauss-Bonnet massive gravity has the advantage of being second 
order in derivatives of metric perturbation, so the Green-Kubo procedure preserves for 
the Gauss-Bonnet massive gravity. 

In the following, after finding a black brane solution in the Gauss-Bonnet massive 
gravity, we calculate the shear viscosity to the entropy density ratio by applying Green- 
Kubo formula. The result shows that this ratio depends on both Gauss-Bonnet coupling 
and the graviton mass and violates KSS conjecture as expected in a higher derivative 
theory. 
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2 Gravity Setup and the Black Brane Solution 


The action of Gauss-Bonnet massive Gravity is as follows Big, 

4 


/= i6 


R + j2+ + m 2 Y Cilii(g, f ) 


L 2 


i=l 


( 6 ) 


C gb = R 2 - 4+ R LWOa R' u ' pa 


pis per J 


where R is the scalar curvature, / is a fixed rank-2 symmetric tensor and m is mass 
parameter. In Eq.©, Cj’s are constants and p are symmetric polynomials of the 
eigenvalues of the 5x5 matrix ICy = \J<f a f av given as mm 


P = [1C] 

U 2 = [1C] 2 - [K 2 ] 

U 3 = [/C] 3 — 3[/C] [/C 2 ] + 2[/C 3 ] 

U 4 = [1C] 4 - 6[/C 2 ][/C] 2 + 8[/C 3 ][/C] + 3 [K 2 ] 2 - 6[/C 4 ] 


The square root in 1C means {y/~A)y{y/~A)\ = Ay and the rectangular brackets denote 
traces. 

We consider the following metric ansatz for a five-dimensional planar AdS black brane, 

ds 2 = ~ T f(r)dt 2 + + r 2 hijdx l dx j , (7) 

L z r z j[r) 

A generalized version of f liu was proposed in BE] with the form = diag( 0,0, c^hij), 
where hij = j?. 

The values of p are calculated as below, 


P 


3c 0 
r ’ 


U 2 


6c| 
r 2 ’ 


P 


6c| 
r 3 ’ 


U 4 = 0 


Inserting this ansatz into the action Eq.© yields, 


I = 


1QttG§ 


d 5 x 


3 N(r) 

~1F~ 


r 4 1 - f(r) + \ gb f{r) + m 2 L z c 0 


°l r 2 ^2 

-y + CQCir + 2c 0 c 3 r 


1 


(8) 

in which 1 denotes derivative with respect to r. Equation of motion is given by variation 
of N(r) [22], 


r 4 ^1 - fir) + X gb f{r) 2 ^j + m 2 L 2 c 0 ^y- + c 0 cir 2 + 2cgC 3 r 


= 0 


f(r) is determined by solving the following equation, 


r A ( 1 - f(r) + A gb f(r) 2 I + ?n 2 L 2 c 0 [ -y + c 0 c 2 r 2 + 2cgC 3 r J = b 4 (9) 
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in which b is an integration constant. This yields two solutions, 


f±( r ) = 


2A, 


gb 


1 ± 


\ 


1 - 4X gb (l - \) - 4m 2 L 2 X gb 


I 2c 0 C 3 C 0 C 2 C 0 Ci 

I o' 




3 r 


( 10 ) 


Since black brane should have an event horizon, we choose /_(r) of the solution above. 
Notice f-(r) is zero at event horizon, r + , so we have, 


7 4 4 

b =r + 


1 + m 2 L 2 


( 2c o c 3 , c o c 2 , cic 0 

' O' 




3 r A 


That’s easy to show that N(r) is constant by variation of f(r) from Eq. 


use the dimensionless variable z = — 

r+ 


~2 2 

ds 2 = - Z ? +^ 


L 2 

with substituting 6 4 in /_(r), 

1 


N f_(z)dt + 


riz 2 3 


L 2 dz 2 


z 2 f~(z) L 2 


5 Z dx t 


i=l 


/-(*) = 


r (-) = 


2A 


gb 


-[l-T(z)] 


( 11 ) 

We 

( 12 ) 

(13) 


\ 


1 + 


4A, 


gb 


1 — z 4 + m 2 L 2 


f 2c[]c3 

V 


(1-Z) + ^P(1-Z2) + ^(1-Z3) 
r± 3r_i_ 


(14) 


In AdS/CFT correspondence, the speed of light in the boundary CFT is simply c = 1. 
So that, for the black brane solution in the asymptotic region we have lini^oo lV 2 /_(z) = 
1 to recover a causal boundary. By applying this criterion we will have, 


n 2 _ 1 + -y/1 4A g b 
2 

The temperature and the Hawking-Bekenstein entropy density follow, 
1 r 1 d 


(15) 


T = 


Nr\ 


2t^ y/g^dr^ 9tt ^ r ~ r + 4vr 
Nr+ 1 1 + m 2 L 2 ( 2cQc3 + 2c(jc 2 + c 0 Cj 


s = 


7 tL 2 
47r 

V 


4 \ 

J d?Xyf^ = 47r (Jj-) > 


3 - 9 ^ + 

n r + r + 


Nr + 
77 L 2 


1 +A 


(16) 

(17) 


where in Eq. (fl6l) . the last equation dehnes parameter A and in Eq. (fT71) . we used 

Teb = 130 b 7 4G - 

We also generalize the metric ansatz to include spherical and hyperbolic, as well as 
planar horizons: 

2 7. r 2 ,. 0 dr 2 


ds z = —(k + —j / {r))N (rf dt 2 + 


k + % fir) 


+ r 2 hijdx l dx J , 


(18) 
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where hijdx l dx J is the following, 

k = +1 : dQ, 2 (metric on S' 3 ), 

1 3 

k = 0 : j2^2( dx i) 2 , 

i =1 

k = —1 : (metric on H 3 ), 

Note that for k = ±1, the above line element has unit curvature. The horizon is 
determined by f(r + ) = — k^-. 


3 Shear Viscosity for the Black Brane 

Let us start with the black brane solution 

o L 2 dz 2 "~ 2 ~ 2 3 


ds 2 = --±JV 2 /df 2 + 

L2 / 




(19) 


Z=1 


in which f = 

r + 

r 2 z 2 ^ 

To find the shear viscosity, we perturb the background metric by /i| = -jjz-tyit, x, z ) 


ds 2 = —±N 2 fdt 2 + 

L2 f 


L 2 dz 2 r+z 2 ( 3 


+ 


L 2 


dx 2 + 2 (f>(t, x, z)dx\dx 2 


( 20 ) 


Z=1 


Plugging in the action ([6]) and introduce Fourier modes of <f>, one finds 

k = (lo, 0,0, q) 


<Kt,x,z)= [ k)e~™ t+ ^ x 


I = 


, (2vr) 2 
— 1 f dzduidq 

V J (2vr) 2 


K(dz4>) 2 - I<24> 2 + d z (K 3 <j) 2 ) 


( 21 ) 


where cj) 2 = 4>(z,k)(j)(z,—k) , 4>(z,—k) = <f>*(z,k ) and it’s functions will be introduced 
in the following. We expand the metric up to second order of 4> and by using the 
general covariance symmetr}Q, the action can be written as follows, 


—1 f dzdudq 


J 


~WF H 


g zz {d z c!)) 2 + g tt (d t <j )) 2 + m 2 ^^-(j) 2 + ... 


( 22 ) 


X A general procedure is given in 
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where, 


K = 

Hg zz = 

II 

£1* 

A 2 

gZZ 

U = 

—Ncqt\ 
4 L 3 

K = 

z 2 f(z- 

H = 

Nr\ 

L 3 T(z) 


El 

L 2 


m 2 U „ 


L 4 


Hg 


it - 


_N 2 f 2 


9 m 2 r 2 ;N 2 U f 


— ° + [3cir + z 2 + 2c 0 c 2 z + 108 cqC 4Z x ] 


(1 - 4A gb )z 6 - X gb L 2 m 2 


( 2cqC 3 2cqC 2 z c 0 ci2 2 

I O ' 


\ 


r + 


(23) 


in which we consider zero momentum q = 0. Notice the surface term K 3 doesn’t 
contribute to the viscosity. Mode equation can be derived from the action m as 
follows, 

K(p" + K'(p' + K 2 cf> = 0 (24) 

Firstly, we solve the mode equation near the horizon. In near horizon limit we have, 


K 2 _ CO 2 

I< ~ B 2 (z- l ) 2 

K' _ 1 

K ~ z - 1 


(25) 

(26) 


By comparing K 2 /K in (1251) and (1231) and ignoring the less singular term proportional 
to m 2 , we obtain, 

5 = !?/'( i) ( 27 ) 

in which ' denotes derivative with respect to 2 . This equation is singular at the horizon 
z = 1. The mode equation at near horizon is, 

<t>"{z) - + Ip (1 - z ) 2 ^ = 0 ’ ( 28 ) 

where its solution is (j)(z ) = (1 — z)P , with 

P = ~ ( 29 ) 

To solve the mode equation (fMj) . we apply the following ansatz, 

cl>(z) = F(zf(l-ph + 0(p 2 )), (30) 


where F(z) = N 2 f(z). Inserting (liffll) in (1241) to first order of (3, we obtain, 

I - Kh" - K'b! + m 2 Uh = 0 . (31) 


KF' 


F 
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It can be easily solved to find, 


(32) 

(33) 

(34) 


(Kyr - Kh')' + m 2 Uh = 0, 

K^r — Kh 1 = A\ — m 2 f~ U(z')h(z / )dz' 
h = log^-A 1 J z ±dz + m 2 f z f Zl U(z 2 )h(z 2 )dz 2 

where A± and A 2 are integration constants. From the above equation, h can be found 
perturbatively. However we need only near horizon and near boundary behavior of h 
as we perform in the Appendix. For m 2 = 0 we have, 

/l( ° ) =l ° g ^~ Al j ~K dZ ( 35 ) 

To be regular at the horizn it is enough to take A\ = K'(1) and this preserves to all 
orders of ?n 2 as shown in the Appendix. We can read the retarded Green’s function by 
expanding the action up to second order of source or the value of field on the boundary, 


S 


1 

2 


d 4 k 

W) 


J a (-k)F a (k-z)J a (k) 



(36) 


then the retarded Green’s function Gn{k ) in momentum space for the boundary field 
dual is given by, 


Gii(k) = -limj a (fc,z)| i, 

£->■0 e 

which can be calculated by the prescription given in 


(37) 


G r (u, 0) = K<j,*d x <f> + K^ 2 = KF-P{ 1 + Ph)PFP(y - h') 


pIa!-eJ U{z')h{z')dz' 


ioj . 

= ~B M 


(38) 


In the first line K^cj) 2 is real and doesn’t contribute to the viscosity, so we drop this term 
in the retarded Green’s function. In the second line the integral term includes terms 
proportional to ci, c 2 and C 4 . A simple power counting indicates that the perturbative 
calculation of h diverges for z —>• 00 unless c\ = c 2 = 0 and C 4 term converges to zero 
as we show in the Appendix. 

Shear viscosity can now be found by using Green-Kubo formula as in the following, 


rj = — lim —QGr(u, 0 ) = = --- 

wow B r + Nf'( 1) 


K\ 1 ) 


We can calculate K'{ 1) as below, 


Hf 1 


K\ 1) = (Hg zz )' | 2= i = (H'g zz + H(g zz )% =1 = H(g zz )'\ z=l = -^ =1 

Then by substituting Ea. (l40l) in Eci. (l39l) . we will have, 

A\ H(z = 1) 


7] = ^ = 


Nr A 


(39) 


(40) 


(41) 
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and the ratio of shear viscosity to entropy density is, 


V 

s 


1 L 3 
47T IVr 3 


H(z= 1). 


( 42 ) 


For the Einstein-Hilbert massive gravity H = y/—g = then the ratio of shear 

viscosity to entropy density is - = j- independent of the graviton mass. 

We apply Eg. (1421) for our black brane, with H given in Eq. (1231) 


! 

s 

1 r 

I^ 1 “ 

4(1 + A)A 9 &] 

(43) 


1 r 

ATrL 2 TX gb 

Nr + J 

(44) 


I^ 1 “ 

where A is in the following, 




A = 

L 2 m 2 f 2cfjc 3 ! 2 cqC 2 ! c 0 ci \ 

(45) 


4 \ r^_ 

r + r + J 

= 

L 2 m 2 CQC3 

2 r 3 


(46) 


where in the last line we put ci = C 2 = 0. For m = 0 we have, ^ [1 — 4A g b] which 

is the well known result for the Gauss-Bonnet gravity. 


4 Conclusion 

In this paper we introduced a new solution to the Gauss-Bonnet massive gravity and 
showed that the lower bound rj/s = 1/47r known as KSS conjecture [26] violates for the 
black brane in this gravity. This is not surprising in the context of higher derivative 
theories m-m- 

Results where exact in both the Gauss-Bonnet coupling and the graviton mass. It 
is interesting that for zero Gauss-Bonnet coupling the graviton mass disappears in the 
viscosity to entropy ratio and the result saturates the KSS bound, 1/47r. On the other 
hand, inclusion of mass to the Guass-Bonnet looks like rescaling the Gauss-Bonnet 
coupling as X gb —>• (1 + A)X g b in the viscosity to the entropy ratio with A a mass 
dependent parameter. Recalling the causal boundary criteria in (jl5l) . it indicates that 
4A g b < 1. Now the rescaling of A g b either decreases or increases the bound on rj/s ratio 
depending on the sign of A. In any case, 1 + A is a positive parameter proportional 
to temperature as given in (fl6l) . 


Acknowledgment Authors would like to thank Ali Imaanpur and Ahmad 
Moradpur for useful discussions. 
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Note added As this article was being completed, we received the preprint [33J 
which has found the same solution as ours. 


Appendix 


Here we study the behavior of h function at the horizon and boundary limits. 
Firstly consider the following perturbative expansion of h in powers of m 2 : 

h(z) = h (0) (z) + m 2 h (1) (z) + m 4 h (2) (*) + ■■■ (47) 

where is given in ()35[) and higher orders can be found by 

= f I " U(z 2 )h^(z 2 )dz 2 (48) 

On the other hand, F and K are regular functions at z = 1, 

KfaK'(l)(z-i). (49) 

Thus to have non-singular h^ at the horizon, A\ is chosen to be A\ = This 

result preserves at higher orders. Notice U is a constant at the horizon, so taking 
A '2 = 1 then h ~ C(z — 1) + • • • with C some constant and 

hW ( z ) ~ / U ( Z ^ C ( Z2 ~ ^ dZ2 = _ 1)2 (50) 

then h^ n \z) = ^ U r\i) ) (£ — l) ri+1 +- ■ ■ • This guarantees that h vanishes smoothly 

near the horizon for any order of m 2 . 

In the same way, we analyze the near boundary behavior of h. For z —>• oo, we have 
K ~ z 3 , F = N 2 f 1, and with A 2 = 1, one finds h 1°) ~ 1 f z 2 . Higher orders of h can 
be derived perturbatively from Eq. (HU). 

Recall U(z) includes terms proportional to ciz 2 , C 2 Z and C 4 Z _1 . It is easy to show 
that positive powers of z in U lead to increasing power of 2 in each step of perturbation 
for h and give a logarithmic divergence in few steps. So we have to take ci = C 2 = 0. 
For C 4 7 ^ 0 the integral converges at the boundary and indeed has zero contribution to 
the viscosity: 


««-7W 

Therefore h^ n \z) r^i l/z 2n+2 and the integral in Eq. (138 p becomes 

f V(z’)h*\z'W ~ f ~ jTj -> 0 


(51) 


(52) 


9 











References 


[1] D. Lovelock, “The Einstein tensor and its generalizations,” J. Math. Phys. 12, 498 
(1971). 

[2] L. A. Gergely, “Brane-world cosmology with black strings,” Phys. Rev. D 74, 
024002 (2006) [hep-th/0603244], 

[3] C. Brans and R. H. Dicke, “Mach’s principle and a relativistic theory of gravita¬ 
tion,” Phys. Rev. 124, 925 (1961). 

[4] M. Akbar and R. G. Cai, “Thermodynamic Behavior of Field Equations for f(R) 
Gravity,” Phys. Lett. B 648, 243 (2007) [gr-qc/0612089]. 

[5] G. Cognola, E. Elizalde, S. Nojiri, S. D. Odintsov, L. Sebastiani and S. Zerbini, “A 
Class of viable modified f(R) gravities describing inflation and the onset of acceler¬ 
ated expansion,” Phys. Rev. D 77, 046009 (2008) [arXiv:0712.4017 [hep-th]]. 

[6] S. H. Hendi, “The Relation between F(R) gravity and Einstein-conformally invari¬ 
ant Maxwell source,” Phys. Lett. B 690, 220 (2010) [arXiv:0907.2520 [gr-qc]]. 

[7] R. G. Cai, Y. P. Hu, Q. Y. Pan and Y. L. Zhang, “Thermodynamics of Black 
Holes in Massive Gravity,” Phys. Rev. D 91, no. 2, 024032 (2015) [arXiv: 1409.2369 
[hep-th]]. 

[8] J. Xu, L. M. Cao and Y. P. Hu, “P-V criticality in the extended phase space of black 
holes in massive gravity,” Phys. Rev. D 91, no. 12, 124033 (2015) [arXiv: 1506.03578 
[gr-qc]]. 

[9] J. M. Maldacena, “The Large N limit of superconformal field theories and super¬ 
gravity,” Int. J. Theor. Phys. 38 (1999) 1113 [Adv. Theor. Math. Phys. 2 (1998) 
231] [hep-th/9711200]. 

[10] O. Aharony, S. S. Gubser, J. M. Maldacena, H. Ooguri and Y. Oz, “Large N field 
theories, string theory and gravity,” Phys. Rept. 323, 183 (2000) [hep-th/9905111]. 

[11] J. Casalderrey-Solana, H. Liu, D. Mateos, K. Rajagopal and U. A. Wiedemann, 
“Gauge/String Duality, Hot QCD and Heavy Ion Collisions,” arXiv:1101.0618 [hep- 
th]. 

[12] Forster D. Hydrodynamic Fluctuations, Broken Symmetry, and Correlation Func¬ 
tions, Reading: Benjamin (1975) 

[13] S. Bhattacharyya, V. E. Hubeny, S. Minwalla and M. Rangamani, “Nonlinear 
Fluid Dynamics from Gravity,” JHEP 0802, 045 (2008) [arXiv:0712.2456 [hep-th]]. 

[14] N. Ambrosetti, J. Charbonneau and S. Weinfurtner, “The Fluid/gravity corre¬ 
spondence: Lectures notes from the 2008 Summer School on Particles, Fields, and 
Strings,” arXiv:0810.2631 [gr-qc]. 

[15] M. Rangamani, “Gravity and Hydrodynamics: Lectures on the fluid-gravity cor¬ 
respondence,” Class. Quant. Grav. 26, 224003 (2009) [arXiv:0905.4352 [hep-th]]. 

[16] L. D. Landau and E. M. Lifshitz, Fluid Mechanics (Course of Theoretical 
Physics,Vol. 6). Butterworth-Heinemann, 1965. 


10 



[17] J. Bhattacharya, S. Bhattacharyya, S. Minwalla and A. Yarom, “A Theory of first 
order dissipative superfluid dynamics,” JHEP 1405, 147 (2014) [arXiv: 1105.3733 
[hep-th]]. 

[18] P. Kovtun, “Lectures on hydrodynamic fluctuations in relativistic theories,” J. 
Phys. A 45 (2012) 473001 [arXiv: 1205.5040 [hep-th]]. 

[19] N. Banerjee and S. Dutta, “Holographic Hydrodynamics: Models and Methods,” 
arXiv: 1112.5345 [hep-th]. 

[20] D.T. Son, Hydrodynamics and gauge/gravity duality, Nuclear Physics B (Proc. 
Suppl.) 192-193 (2009) 113-118. 

[21] Joseph I. Kapusta , Charles Gales; “finite- temperature held Theory Principle and 
applications” Cambridge university press 2006. 

[22] R. C. Myers and B. Robinson, “Black Holes in Quasi-topological Gravity,” JHEP 
1008, 067 (2010) [arXiv:1003.5357 [gr-qc]]. 

[23] R. G. Cai, Z. Y. Nie and Y. W. Sun, “Shear Viscosity from Effective Couplings 
of Gravitons,” Phys. Rev. D 78, 126007 (2008) [arXiv:0811.1665 [hep-th]]. 

[24] G. Policastro, D. T. Son and A. O. Starinets,“From AdS/CFT correspondence to 
hydrodynamics,” JHEP 0209, 043 (2002) [hep-th/0205052]. 

[25] P. Kovtun, D. T. Son and A. O. Starinets, “Holography and hydrodynamics: 
Diffusion on stretched horizons,” JHEP 0310, 064 (2003) [hep-th/0309213]. 

[26] P. Kovtun, D. T. Son and A. O. Starinets, “Viscosity in strongly interacting 
quantum held theories from black hole physics,” Phys. Rev. Lett. 94, 111601 (2005) 
[hep-th/0405231]. 

[27] V. Jahnke, A. S. Misobuchi and D. Trancanelli, “The Chern-Simons diffusion rate 
from higher curvature gravity,” Phys. Rev. D 89, 107901 (2014) [arXiv: 1403.2681 
[hep-th]]. 

[28] M. Brigante, H. Liu, R. C. Myers, S. Shenker and S. Yaida, “Viscosity 
Bound Violation in Higher Derivative Gravity,” Phys. Rev. D 77, 126006 (2008) 
[arXiv:0712.0805 [hep-th]]. 

[29] M. Brigante, H. Liu, R. C. Myers, S. Shenker and S. Yaida, “The Viscosity Bound 
and Causality Violation,” Phys. Rev. Lett. 100, 191601 (2008) [arXiv:0802.3318 
[hep-th]]. 

[30] I. P. Neupane and N. Dadhich, “Entropy Bound and Causality Violation in Higher 
Curvature Gravity,” Class. Quant. Grav. 26, 015013 (2009) [arXiv:0808.1919 [hep- 
th]]. 

[31] A. Bhattacharyya and D. Roychowdhury, “Viscosity bound for anisotropic su- 
perhuids in higher derivative gravity,” JHEP 1503, 063 (2015) [arXiv: 1410.3222 
[hep-th]]. 

[32] D. Roychowdhury, “Hall viscosity to entropy ratio in higher derivative theories,” 
JHEP 1410, 15 (2014) [arXiv: 1408.5095 [hep-th]]. 

[33] S. H. Hendi, S. Panahiyan and B. E. Panah, “Charged Black Hole Solutions in 
Gauss-Bonnet-Massive Gravity,” arXiv: 1507.06563 [hep-th]. 


11 



